DIGITAL MEASUREMENTS OF SPREAD SPECTRUM CLOCKING 



Background Of The Invention 

The FCC in the United States and other regulatory bodies in foreign jurisdictions have become 
more aggressive in setting performance standards for the radiation of EMI (Electro-Magnetic 
Interference). Digital equipment that operates according to clock signals, such as computers and their 
mass storage peripherals, have been identified as a source of EMI. Historically, a virtue of a clock signal 
was frequency stability; many data transmission systems have crystal controlled clocks. The combination 
of this frequency stability and the fast edges of the clock and its associated signals combine to produce 
pronounced collections of harmonics that can be the cause of EMI. For various reasons, it might not be 
technically or economically feasible to prevent the EMI through shielding. One of the techniques that 
has been developed to reduce EMI from such systems is to apply spread spectrum techniques to the clock 
system. SSC (Spread Spectrum Clocking) is a technique where the frequency of the clock signal is made 
to vary with time in some controlled manner. 

For example, there is a Serial ATA standard for hard disc drives that specifies that the clock 
frequency of 75 MHZ (or 66 MHZ) vary with respect to time downwards and then back up again by one 
half percent, and at a rate that varies from 30 KHz to 33 KHz. That is, the clock frequency undergoes 
a varying degree of Frequency Modulation (FM). The idea is to distribute the harmonic content 
throughout a portion of the spectrum to reduce its power level at any particular frequency. 

The FCC and the other regulatory bodies have certification requirements that require the 
manufacturer to measure the performance of the system to ensure that it meets the regulatory 
requirements. Items of interest in a system using SSC include the degree of power reduction at certain 
frequencies (e.g., the nominal clock frequency without FM) and a description of the modulating FM 
signal itself. Now it is not the case that there is no test equipment available to make such measurements; 
there is. But the circuitry that needs to be tested often resides on an IC (Integrated Circuit) die that will 
be tested on a system that is designed to probe, power up, stimulate and test die locations on undiced 
wafers. They are well suited to the task of ordinary functional testing of digital operation, but do not 
incorporate the kinds of capabilities found in the type of laboratory grade bench test equipment (flexible 
counters, fast digital oscilloscopes and full featured spectrum analyzers) that would ordinarily be used 



to make measurements related to SSC. What is more, it would be electrically troublesome, not to 
mention expensive, to connect such bench test equipment to a wafer or die. Yet it is most desirable in 
high volume production situations to test parts as early as possible, so as to avoid the expense of building 
assemblies with bad parts, since the cost of finding and replacing a bad part goes up rapidly as assembly 
5 fabrication proceeds. The parts we are interested in do get otherwise tested by a digital IC tester such as 

the Agilent 93000 SOC System (SOC stands for System On a Chip). However, depending upon how any 
particular such system is configured, the tester might not be equipped with the hardware needed to 
perform measurements of RF spectrum and FM characteristics. The Agilent 93000 SOC System can 
indeed be equipped with such measurement hardware, but to do so involves a certain expense. Such an 
10 expense might be viewed as a lamentable inefficiency if it were incurred only for the measurement of 

parameters related to SSC. It would desirable if there were a way that the measurement capabilities of 
the Agilent 93000 SOC System, and of other similar systems, could be easily and at a low cost, extended 
to measure the spectrum, deviation and rates of deviation for an SSC operated system. What to do? 

Summary Of The Invention 

15 A solution to the problem of measuring the RF spectrum, and FM characteristics of the clock 

signal in an SSC operated system is to take advantage of one of the one or more high speed single-bit 
test channels that can be present in the Agilent 93000 SOC System. Such a channel can be configured 
to measure, at up to 2.5 Gbps, the logical value of an input signal applied thereto. The SSC signal of 
interest is applied to one of the high-speed single-bit test channels. A conventional FFT (Fast Fourier 

20 Transform) is performed on the captured data to discover the aggregate nature of the distributed spectral 

components, which can then be compared with associated specifications. The captured data is applied 
to another algorithm to find the FM modulation profile. That algorithm involves operating on captured 
digital data to perform a frequency translation operation (equivalent to heterodyning), a first filtering 
operation, a discrete differentiation operation and raw phase extraction that includes a summation of two 

25 components and a root square operation, finally followed by a second filtering operation to produce an 

FM modulation profile. The algorithm for finding the SSC modulation profile does not require the high- 
speed digital channel to meet Nyquist sampling requirements. 
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Brief Description Of The Drawings 

Figure 1 is a simplified flow chart of how SSC measurements are made in accordance with an 
aspect of the invention; 

Figures 2A and 2B are spectrum diagrams useful in appreciating the operation of the simplified 
5 flow chart of Figure 1 ; and 

Figures 3 A and 3B are respective phase and frequency diagrams illustrating the recovery of an 
FM modulation profile from an SSC signal that has been processed to extract phase changes as a 
function of time. 

Description Of A Preferred Embodiment 

1 0 Refer now to Figure 1 , wherein is shown a simplified flow chart 1 that includes an algorithm for 

obtaining the FM parameters of an SSC signal from captured digital data that does not necessarily meet 
Nyquist criteria. The first step 2 in the flow chart 1 is the capture of the SSC data using the "digital 
channel" capability of the Agilent 93000 SOC System. In this type of measurement an SSC signal is 
applied to a comparator having a threshold and the logical results of TRUE and FALSE appear as a 

1 5 stream of digital data at that matches the SSC signal. A digital 11 1 " will represent TRUE and a digital "0" 

will represent a FALSE. The threshold for the comparator is set at about the mid-point of the total 
expected excursion of the SSC signal. At some regular and pre-determined rate the digital values in this 
single bit stream are captured in a array in memory. The rate of capture will not be synchronous with the 
SSC signal, but will instead be at a regular and steady rate determined by the sampling mechanism in 

20 the digital channel of the tester. Logic State Analyzer and Timing Analyzer hardware come to mind as 
representative of the sort of data capture hardware being described here. Now, if the data capture were 
synchronous with the SSC signal, then the captured data would lose its time reference and also its SSC 
(FM) character. However, a steady rate of capture that is at least as fast as the fastest clock rate (which 
we shall call co 0 ) is guaranteed to preserve the underlying FM properties, even though the captured data 

25 may from time to time be "missing" a bit. That is, the regular nature of sampling may, at rates lower than 
twice the highest SSC signal rate, occasionally fail to register a transition. The result is that a string of 
consecutive Ts might be missing a 0, or a string of consecutive 0's might be missing a 1 . The algorithm 
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that we will disclose for recovering the FM modulation profile tolerates this minor wart, and produces 
the correct answer, anyway. These measured Ts and O's are held in a first measurement array. 

It will further be appreciated that, since the algorithm begins with a threshold operation, the 
signal whose FM profile is to be measured might be a sinusoidal FM signal having equal plus and minus 
5 peak excursions about ground, rather than a digital signal referenced to ground. In such a sinusoidal case, 
the threshold would be, say, greater than or equal to zero volts. Furthermore, the modulating signal (i.e., 
the signal whose modulation profile is of interest) need not be sinusoidal; it might, for example, be 
triangular, or even something else, such a as a saw tooth. Also, and in either case, it is desirable to 
capture at least three to five cycles of the SSC signal, to allow for settling of certain filtering steps to be 

1 0 described in due course. 

To continue, then, if base band and harmonic spectrum information are desired, then the captured 
data is (at 3) operated upon with a suitable (and conventional) DFT (Discrete Fourier Transform) 
process, whereupon the results can be analyzed and classified (at 4), after which the spectrum data can 
be examined and compared to specified limits (5). This right-hand fork of the flow chart 1 may require 

1 5 that the capture rate of the digital channel meet certain Nyquist criteria. 

Our main interest is in the left-hand fork of the flow chart 1. Here there are shown a series of 
steps 6-10 that operate on captured SSC data acquired back at the first step (2). The captured data might 
be the same data that was used by the right-hand branch of the flow chart, or it might be a data from a 
separate instance of data acquisition to supply data in the first measurement array. 

20 At this point the reader should refer simultaneously to Figures 1 , 2A and 2B. As shown in Figure 

2A, the SSC signal under consideration has a certain spectrum 1 1 (we have not attempted to indicate any 
harmonics - just the base band, or fundamental, spectrum itself). It is the FM modulation profile of that 
spectrum 1 1 that we wish to measure. The TRANSLATE SPECTRUM step (6) of the flow chart 1 is 
a software implementation of a heterodyne operation that produces upper and lower side bands (12 and 

25 14, respectively in Figure 2B) that are not only displaced from ?b 0 , but that are also reflected (or reversed 

in "direction 1 '). In order to perform this computation we employ the expedient of taking a logical 1 to 
represent the number one, and taking a logical 0 to represent a the number zero. It will be appreciated 
that the modulation profile of the sidebands is, save for the reflection, the same as for the original 
spectrum 1 1. The TRANSLATE SPECTRUM step (6) operates on the first measurement array and 
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produces in a second translation array a (numerical) string of plus or minus ones and zeros that 
represents the superposition of the two side bands. 

The next step (7) in the algorithm is the application to the data in the second translation array 
of a FIRST FIR (Finite Impulse Response) FILTER. Its purpose is to act as a low pass filter that removes 
5 the upper side band. The output of the FIRST FIR FILTER is a third array that would normally have both 

original sine and original cosine sections that are populated with fractional values for each sample, save 
that the values of the cosine terms in the cosine section will be all zeros, and need not actually exist as 
memory having those zero values. (The absence of any non-zero cosine terms arises from setting the 
sampling threshold at the mid-point of the excursion of the SSC signal, and at a rate we shall call co 0 in 

10 Figure 2 A and 2B. The computations of the FIRST FIR FILTER may be simplified to reflect this fact. 

If the sampling were accomplished at twice that rate then there would be two interleaved original 
sequences of sine and non-zero cosine data. That is, if sampling were done at a rate of 2co 0 , then the 
result would be construed as two interleaved sections that were each sampled at co 0 and that are separated 
in phase by 90°. What is of note is that the algorithm shown below still works if the data is sampled at 

15 justco 0 .) 

Following the FIRST FIR FILTER step (7) is a DISCRETE DIFFERENTIATION AND RAW 
PHASE EXTRACTION STEP (8). The discrete differentiation takes the differences between successive 
original sine terms in the third array to operationally discover the rate of change of those original sine 
terms, which is, as is well known, cosine information. If the original cosine terms are not zero, then the 

20 same discrete differentiation is also performed for the original cosine terms to find derived sine terms, 

although this step maybe omitted if the original cosine terms are zeros. The original sine terms and their 
derived cosine terms, and the original (and assumed-to-be-zero) cosine terms and their derived sine 
terms can all be combined to find phase information, as shown below. The change in the phase (Acp(f)) 
can then be used to further find the change in frequency (Acg(0), also as shown below. It is, of course, 

25 the change in frequency that describes the FM modulation profile. 
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The next step is a SECOND FIR FILTER (9). This accepts a fourth array of consecutive Aco(7) 
values from step (8) and smooths it to become a fifth array of consecutive values that the FM modulation 
profile. An example of what such data is like is shown in the graph 15 of Figure 3 A. The filtered result 
is not necessarily perfect, a typical graph 16 of which is shown in Figure 3B; location 17 is a minor 
imperfection that arises from nature of the binary sampling technique, and might be ameliorated by 
sampling at a higher rate or by more stringent filtering. Note also that the SECOND FIR FILTER needs 
some time to settle (i.e., note that the extreme left-hand edge of the graph is misshapen). This is a 
principle reason why the CAPTURE step (2) ought to acquire three to five cycles of the SSC signal. It 
will be appreciated that the recovery of different modulation profiles for the SSC signal might benefit 
from corresponding adjustments to, or changes to the nature of, the SECOND FIR FILTER. 

Subsequent to the SECOND FIR FILTER step the sought-after FM modulation characteristics 
are available and can be examined (10). That is, the fifth array can be treated as a potential graph (e.g., 
1 6 in Figure 3B) that is in tabular form, and inspected for maximum and minimum ordinates (which after 
suitable scaling represents frequency values) along the abscissa (time). 



The following APPENDIX is pseudo code (essentially a complete C++ program) that implements 
the steps of the left-hand fork of the flow chart of Figure 1. In one embodiment this program executes 
on the Agilent 93000 SOC System. It is well commented and is believed to be, in light of the foregoing 
Specification, self-explanatory to those who are skilled in the C++ programming language and also in 
5 the basics of digital signal processing techniques, such as the Discrete Fourier Transform. The line 
numbers that appear to the left of the listing are not part of the program and have appended merely to 
facilitate references in this Specification to locations in the text of the APPENDIX. 
In particular: 

Lines 20 - 56 correspond to step (2) of the flow char of Figure 1; 
10 Lines 57-72 correspond to step (6) of the flow char of Figure 1; 

Lines 73 - 93 correspond to step (7) of the flow char of Figure 1; 

Lines 94- 123 correspond to step (8) of the flow char of Figure 1; 

Lines 124- 137 correspond to step (9) of the flow char of Figure 1; 

Lines 138- 165 correspond to step (10) of the flow char of Figure 1; and 
15 Lines 166 - 220 are a subroutine for the FIR filters of steps (7) and (9). 
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APPENDIX 



1 // 

2 //- DESCRIPTION: 

3 II- 

4 //- This pseudo code (which is almost a complete C++ program) shows 

5 //- the extraction of SSC frequency modulation information based on 

6 //- captured single bit digital data which might look like, for 

7 //- example: 101010101010001010101010110110101010101010 . 

8 II- 

9 II- This example routine executes on an Agilent 93000SOC ATE 

10 //- system. It captures data on a single pin under consideration, 

11 //- which will be the SSC output pin or transmitter pin which 

12 //- outputsan SSC signal at high speed. It uses a standard 

13 //- subroutine call available in Agi lent93000SOC SW, through which 

14 //- captured data will be retrieved into the workstation. 

15 //- Subsequently, this algorithm will begin to operate on this 

16 //- captured single bit digital data to extract frequency 

17 //- modulation parameters. 

18 II- 

19 // 

20 // 

21 //- A) Digital capture test: 

22 //- 

23 //- A.l) 

24 //- First an array variable Captured_Data of type ARRAY_D 

25 //- is defined. ARRAY_D ' s an one-dimensional array data of 

26 //- DOUBLE type. The system automatically allocates memory 

27 //- by the number of elements assigned to "Captured_Data" . 

28 //- The number of array elements is checked and is automatically 

29 //- resized if necessary. 

30 II- 

31 //- A. 2) 

32 //- DIGITAL_CAPTURE_TE ST ( ) is a standard API ( Application 

33 //- Program Interface) call of 93000 SOC SW to run the test in the 

34 //- digital capture mode, and capture data on the specified 

35 //- digital pin for SSC clocking. The captured data is saved in 

36 //- a memory associated with the pin ("tester pin channel memory'). 

37 II- 

38 //- A. 3) 
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39 //- Retrieve data from the vector memory 

40 //- VECTOR("VectorVariable") .getVectors(NumberSamples) is 

41 //- another standard API call for retrieving data from the 

42 //- tester pin channel memory back into the array variable 

43 //- "Capture_Data" defined before. 

44 //- "VectorVariable" is data in tester pin channel memory 

45 //- associated with the SSC signal pin, and 

46 //- "NumberSamples" is the number of samples specified. 

47 // 

48 ARRAR_D Capture_Data ; 

49 DIGITAL_CAPTURE_TEST() ; 

50 Capture_Data=VECTOR(VectorVariable) .getVectors(NumberSamples) ; 

51 // 

52 //- A. 4) 

53 //- Put the captured data to graphic window for checking. 

54 //- PUT_DEBUG is an standard API to perform this task. 

55 // 

56 PUT_DEBUG("Pin" , "RawData" ,Capture_Data) ; 

57 // - 

58 //- B) Spectrum Translation with sin: 

59 II- 

60 //- Only sin data will be calculated, since all cos data is 

61 //- is equal to zero. 

62 // 

63 ARRAY_D AutoWSi n , AutoWCos ; 

64 cap_size=Capture_Data . size() ; 

65 AutoWSi n. resize (cap_size) ; 

66 AutoWCos. resize(cap_size) ; 

67 INT loopcount , cap_si ze ; 

68 for (loopcount=0;loopcount<cap_size;loopcount+=2) 

69 { 

70 AutoWSi n [loopcount ]=Capture_Data[ loopcount] 12 ; 

71 AutoWSi n [loopcount+1] =-Capture_Data [loopcount+1] 12 ; 

72 } 



73 // 

74 //- C) First FIR Filtering: 
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75 II- 

76 //- C.l) 

77 //- First FIR generation: 

78 //- dsp_LPF_FIR is a subroutine for creating coefficients 

79 //- of a finite impulse response filter. 

80 //- 

81 //- C.2) 

82 //- Applying this filter to data: 

83 //- DSP_CONVOL is a standard API call to perform convolution 

84 //- operation of two arrays and. result is saved in FIRWSin 

85 //- and FIRWCos. FIRWCos in fact are all zero. 

86 // 

87 INT Nfir=1024; 

88 DOUBLE dLPBW=0.05; 

89 ARRAY_D dFIR; 

90 dsp_LPF_FIR(Nf i r .dLPBW.dFIR) ; 

91 ARRAY_D FIRWSin , FIRWCos ; 

92 DSP_C0NV0L (dFIR, AutoWSi n .FIRWSin) ; 

93 DSP_C0NV0L (dFIR, AutoWCos, FIRWCos) ; 



94 // - 

95 //- D) Discrete Differentiation and Raw Phase Extraction 

96 //- 

97 //- D.l) Discrete Differentiation based on FIRWSin(i)-FIRWS(i-l) 

98 //- and FIRWCos (i )- FIRWCos (i -1) where FIRWSin and FIRWCos 

99 //- are the sinusoidal and cosinusoidal components of the 

100 //- data array from the first FIR filter result. 

101 II- 

102 //- D.2) Squaring of each differentiation result 

103 II- 

104 //- D.3) Summation of these two squaring results and root 

105 //- squaring afterwards 

106 // 

107 ARRAY_D SSC_Value (Capture_Data . si ze () ) ; 

108 DOUBLE Tmp_Value; 

109 for (count=0;count<points;count++) 

110 { 

111 if (count==0) 

112 { 

113 SSC_Value[count]=0.0; 
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114 } 

115 else 

116 { 

117 Tmp_Value=(FIRWSin [count] -FIRWSin[count-l] )* (FIRWSin [count] -FIRWSin 

118 [count-1]); 

119 Tmp_Value+=(FIRWCos [count] -FIRWCos [count-1] )* (FIRWCos [count] -FIRWCo 

120 s[count-l]); 

121 SSC_Value [count] =sqrt(Tmp_Value)/Step/2/M_PI*2; 

122 } 

123 } 

124 // - 

125 //- E) Second FIR Filter: 

126 II- 

127 II- E.l) 

128 //-Second FIR Generation: 

129 //- 

130 //- E.2) 

131 //- Filtering with the second FIR filter. 

132 //- After the filtering the desired FM parameters are present. 

133 // - 



134 ARRAY_D SSC_Demod; 

135 dLPBW=0. 000005; 

136 dsp_LPF_FIR(Nf*i r .dLPBW, dFIR) ; 

137 DSP_CONVOL(dFIR,SSC_Value,SSC_Demod) ; 



138 // 

139 //- F) Analysis of the result 

140 II- 

141 //- Find the frequency data from the changing phase data. Find 

142 //- the minimum and maximum values of the frequency modulation. 

143 // 

144 // 

145 //- DSP_MINMAX is a standard API call of 93000 SOC SW 

146 //- to find the min and max value in an array. The variables 

147 //- min and max are proportional to the min and max frequency. 

148 // - 

149 DOUBLE min, max; 

150 DSP_MINMAX(SSC_Demod,&min,&max) ; 



- 11 - 



152 //- DSP_EDGE is a standard API call of 93000 SOC SW 

153 //- to find the element location of rise/fall in an array 

154 //- with specified threshold. 

155 //- First find edgel at 50% threshold and then find the next 

156 //- edge2 at 50% threshold. The difference is then between 

157 //- edge2 and edgel and result is saved in period. 

158 // - 

159 DOUBLE edgel, edge2; 

160 DOUBLE period; 

161 int start; 

162 start=1000; nTet . s 

163 DSP_EDGE(SSC_Demod, &edgel, start, 0.5*(max-min)+mm; RISE); 

164 DSP_EDGE(SSC_Demod, &edge2, edgel, 0.5*(max-min)+min; RISE); 

165 period=edge2-edgel; 

166 // 

167 // End of the main routine. 

168 // "'" 

169 // " 

170 // Subroutine for finite impulse response filter coefficients: 

171 // 

172 void dsp_LPF_FIR( 

173 INT Nfir, /* Length of FIR ' 

174 DOUBLE dBW, /* Pass band width of LPF (1.0 = Fdgt) */ 

175 ARRAY_D & dFIR /* LPF FIR array data */ 

176 ) 

177 { 

178 INT i .Nsp.Nc.Ncycles; 

179 ARRAY_COMPLEX CSp.CWave; 

180 ARRAY_D dWavel , dWave2 ; 

181 ARRAYJ) dSpl,dSp2; 

182 DOUBLE dX.dNc.dP; 

183 dNc=Nf ir*dBW; /* Passband Ratio (relative to dFdgt) / 

184 Nc=(int) (dNc+0.5) ; 

185 Nsp=Nfir/2; 

186 CSp. resize(Nfir) ; 

187 for (i=0;i<=Nc;i++) 

188 { 

189 CSp[i] .real()=1.0; 
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190 CSp[i] . imag()=0.0; 

191 } 

192 for (i=Nc+l; i<=Nsp; i++) 

193 { 

194 CSp[i] .real()=0.0; 

195 CSp[i] .imag()=0.0; 

196 } 

197 for (i=l;i<Nsp;i++) 

198 { 

199 CSp[Nfir-i] .real()= CSp [i ] . real () ; 

200 CSp[Nfir-i] . imag()=-CSp[i] .imag() ; 

201 } 

202 DSP_IFFT(CSp,CWave) ; 

203 dFIR.resize(Nfir) ; 

204 for (i=0;i<Nfir;i++) dFIR[i ] =0. 5*CWave [i ] . real () ; 

205 for (i=0;i<Nsp;i++) 

206 { 

207 dX=dFIR[i] ; 

208 dFIR[i]=dFIR[Nsp+i] ; 

209 dFIR[Nsp+i]=dX; 

210 } 

211 /* Hanning weighting onto the FIR */ 

212 for (i=0;i<Nfir;i++) 

213 { 

214 dX=0. 5* (1.0-cos (2. 0*M_PI* (double) i/ (double) (Nfir-1))) ; 

215 dFIR[i]=dFIR[i] *dX; 

216 } 

217 /* Scaling */ 

218 dX=0.0; for (i=0 ; i <Nf i r ; i++) dX=dX+dFIR[i ] ; 

219 DSP_MUL_SCL(1.0/dX,dFIR,dFIR) ; 



220 } 
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